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We discuss how electromagnetic properties provide useful tests of the nature of resonances, and 
we study these properties for the A r *(1535) which appears dynamically generated from the strong 
interaction of mesons and baryons. Within this coupled channel chiral unitary approach, we evaluate 
the Ax ii and S1/2 helicity amplitudes as a function of Q 2 for the electromagnetic A r *(1535) — > j*N 
transition. Within the same formalism we evaluate the cross section for the reactions jN — > r)N. 
We find a fair agreement for the absolute values of the transition amplitudes, as well as for the 
Q 2 dependence of the amplitudes, within theoretical and experimental uncertainties discussed in 
the paper. The ratios obtained between the S1/2 and Ai/ 2 for the neutron or proton states of the 
-/V*(1535) are in qualitative agreement with experiment and there is agreement on the signs. The 
same occurs for the ratio of cross sections for the rj photoproduction on neutron and proton targets 
in the vicinity of the AT* (1535) energy. The global results support the idea of this resonance as 
being dynamically generated, hence, largely built up from meson baryon components. However, the 
details of the model indicate that an admixture with a genuine quark state is also demanded that 
could help obtain a better agreement with experimental data. 

PACS numbers: 14.20.Gk,13.40.Gp,12.39.Fe 
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I. INTRODUCTION 

The traditional picture of baryons as being made from 
three constituent quarks EI is giving rise, in some cases, 
to more complicated structures. One of the ideas which 
has gained strength in recent times is that low lying reso- 
nances of J p = l/2~, 3/2~ seem to be well represented in 
terms of states which are generated by the meson baryon 
interaction in L = 0; in the 1/2 - case from the inter- 
action of the octet of mesons of the it with the octet of 
baryons of the p @, i, i, i, i, 0, 0, H and in the 3/2" from 
the interaction of the same mesons with the decuplet of 
baryons of the A(1232) EI E2- The A(1405), which 
actually comes as two poles in chiral theories Q, with 
this two-pole structure supported by experiment [l2l ]. 
has been for long thought of as a kind of meson baryon 
molecule of the KN and 7rE states EH E3 1 a structure 
similar to that provided by the chiral approaches men- 
tioned above. The A/* (1535) is one more resonance that 
appears in the two octets and one singlet of dynamically 
generated resonances coming from the interaction of the 
octet of mesons of the n with the octet of baryons of 
the p Q. In fact, it was noted earlier in [l5| . before 
the systematics of [f| was established, that the inter- 
action provided by chiral Lagrangians put as kernel of 
the Lippmann Schwingcr equation generated this reso- 
nance, which also appears in other work along similar 
lines PIE!- 

The AT* (1535) plays an important role in all processes 
of 77 production since it couples very strongly to rjN . This 
feature is actually provided automatically by the chiral 
theories, one of the points of support for the nature of 
this resonance as being dynamically generated. A recent 
study of the model dependence of the properties of this 



resonance is seen in [181 ] . 

From the point of view of a dynamically generated 
resonance the AT* (1535) leads to fair descriptions of the 
ttN — > rjN and 7 A" — > i]N reactions EH EtI Il9j] and pro- 
duces reasonable numbers for the rjN scattering lengths 
@,E3- Yet, it has been argued that one of the important 
tests of the nature of a resonance is its electromagnetic 
form factors. Indeed, a meson baryon resonance should 
get the Q 2 dependence basically from the meson cloud 
(we take as usual Q 2 = —k 2 with k being the photon 
momentum). If this is a pion, this light particle has a 
fairly large extent in the wave function, as a consequence 
of which, the form factor of the resonance should fall rel- 
atively fast compared to ordinary quark models which 
confine the quarks at smaller distances. This is also the 
case for the proton at small Q 2 , due to its meson cloud, 
which stabilizes later on at larger values of Q 2 where the 
quark c omp onents take over, as shown in chiral quark 
models [2(| [2l|, H2J. We shall see that something spe- 
cial happens for the AT* (1535), but in any case this is a 
very stringent test, since the chiral theory provides the 
normalization and the Q 2 dependence for the different 
transition form factors without any free parameter, once 
the parameters used in nN scattering with its coupled 
channels are fixed to scattering data. 

Radiative decays of resonances from the point of view 
of their dynamically generated nature have been ad- 
dressed in [13 for the A(1520), in [H for the A(1700) 
and in [25| for the two A(1405) states. It concerns the 
decay of the resonances into a baryon and a real photon. 
Some work with virtual photons from this point of view is 
done in 0, [26[ for the electroproduction of 77 in the vicin- 
ity of the A"* (1535) resonance. Meanwhile, experimental 
analyses have succeeded in extracting the helicity transi- 
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tion form factors for A r *(1535) — > Nj A1/2 and S1/2, for 
both N = p,n, in a relatively wide range of Q 2 values 

We evaluate these form factors from the point of view 
of the A r *(1535) as a dynamically generated resonance. 
For that purpose we shall extend the formalism of [13, [2f| 
to virtual photons. The new formalism requires changes 
from the real photon case, but it is rewarding since it 
provides much more information, replacing the helicity 
transition amplitudes by functions of Q 2 and adding the 
new 51/2 transition form factor which only plays a role 
for virtual photons. Hence, there is far more information 
to test the predictions of the model. 

From the quark model point of view there has also 
been much work done on these helicity form factors 
[H, HI M, HI . A comparison of their 
prediction with experiment plus a compilation of results 
from different experiments can be seen in 0, H3, 0]- 
There are appreciable differences from one quark model 
to another and relativistic effects seem to be important, 
particularly in the S1/2 helicity transition form factor. 
It should be noted that some of the models, particularly 
those incorporating relativistic effects 0, 0, pro- 
duce a fair agreement with data, in particular a good 
description of the Q 2 dependence of the form factor. 

In our approach, the quarks enter through the meson 
and baryon components of the resonance and the Q 2 de- 
pendence is tied to the meson and baryon form factors, 
which we take from experiment, plus the particular Q 2 
dependence of the loop functions from the meson baryon 
coupled channels that build up the resonance. Thus, the 
final Q 2 dependence is a non-trivial consequence of chiral 
dynamics, which provides the coupling of the resonance 
to open and closed channels, the Q 2 dependence of the 
different loops, and the form factors of the mesons and 
baryons, particularly the mesons, as we shall see. 

The results that we obtain are in fair agreement with 
experiment for both charged states and for the two tran- 
sition form factors, hence providing extra support for the 
nature of the iV*(1535) as being largely made from the 
interaction of meson baryon coupled channels. 

Some deficiency in the Q 2 dependence at large Q 2 
could be an indication of a mixture of the meson baryon 
components with a genuine quark component, which is 
also indicated by particular details of the chiral approach 
that we shall mention below. 



II. HELICITY AMPLITUDE 

We consider the production reaction of the 7V(1535) 
resonance (J p = l/2~) by in-elastic electron- nucleon 
scattering as shown in FigQ] The 7V(1535) is created by 
exchange of a virtual photon carrying momentum k. The 
initial N and final N* momenta and masses are denoted 
by (pi, Mm) and (P, Mm*), respectively. The energy mo- 




FIG. 1: Kinematics of the electroproduction of iV(1535). 

mentum conservation reads 

P = n + k. (1) 

There are two independent amplitudes for the electro- 
transition from J p = l/2 + to I /2~, A1/2 and 5*1/2, which 
are defined in terms of the transition electric current J M 
by 

A 1/2 = ] ^- e (N*,J z = ±\e™J»\N,S, = -\) (2) 

with the fine structure constant a = e 2 /4ir, the en- 
ergy equivalent to that of a real photon qn = (W 2 — 
M N )/{2W) and the photon-nucleon center-of-mass en- 
ergy W = VP 2 . The polarization vectors of the photon, 
e^, are given by 

4 = -L(0,Tl,-i,0) (4) 
e £ = -L=(fc,0,0,-fc o ) (5) 

VQ 

with Q 2 = — k 2 , where we take the CM momenta k and 
pi along the z axis. 

Let us discuss the general expression of the transition 
current J M in the relativistic formulation. First of all, we 
recall the equation of motion for the initial nucleon 

(jfc - M N ) Ui (pi) = (6) 

where Ui(pi) is the Dirac spinor for the initial nucleon 
normalized by 

j Ej + M N ( 1 \ 

For the final N* , we assume the pole dominance, so that 
we again have 

(f-M N ,)u f (P)=0 (8) 

where Uf(P) is the N* Dirac spinor and Mm- denotes 
the real part of the N* mass. In the calculations of the 
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hclicity amplitude, the Mjv* is chosen as the N* energy 
in the final state, W ~ 1535 MeV. 

It follows that the terms involving the 7 matrix in J ^ 
are only of the form 7-6, since we can move fa and p 
through 7^7^ — Ivln = 9nu to either the left or right 
end in the amplitudes, and they can be replaced by the 
masses by means of Eqs. © and ([5]). The term involving 
ft can also be replaced by the momentum conservation 
ft = p — fa . Thus, Lorentz invariance and momentum 
conservation |T]) require the transition current J M to be 
written, in general, by the following three Lorentz scalar 
amplitudes: 



J" = (M^ + M 2 P" + M 3 k») l5 



(9) 



The gauge invariance k • J = 0, tells us that there 
are only two independent amplitudes among these three 
amplitudes, Mi, giving the following relation: 



(M N , + M N )Mi + k ■ PM 2 + k 2 M 3 = 



(10) 



Using the transition current ©, we evaluate the helic- 
ity amplitudes, A\j 2 and S1/2, in the rest frame of the 
iV(1535) resonance. After some algebra, the helicity am- 
plitudes are written in terms of the amplitude M 2 and 
M 3 by 
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(11) 



\k\ 



2M N e M N * + M N 
x {M N ,M 2 + (M N , - M N )M 3 ) (12) 

The transition current Q can be written equivalcntly 
in the CM frame as 



Ei + M 



N 



2M 
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Mi 



(Eil 
M 3 



+ 



M 2 



M N )W E,+M_ 



N 



P»o ■ k 



-k' 



E t + M N 
TVl^V 1 + M^P^a ■ k 



M^ R k"a ■ k 



(13) 
(14) 



where cr M = (0, a) and we take the CM frame P M = 
(W, 0). Then the helicity amplitudes are written in terms 
of the amplitudes defined above, Mf R 



as 



A1/2 = 
S\/ 2 = 



f^LV2(k-PMr 



2ira 1 

qR e 
2-ixa -|fc| 

1R e 



k 2 M* R : 



(WM 



NR. 
2 



k°M^ R ) 



(15) 
(16) 



with the gauge invariance condition for the nonrelativis- 
tic amplitudes 



+ M^k ■ P 



NR 1 



M 3 R k 2 



0. 



(17) 



III. EVALUATION OF THE TRANSITION 
FORM FACTORS 

A. Model of iV(1535) and photon coupling 

In our approach, the A^(1535) resonance is dynami- 
cally generated in the s-wave meson baryon scattering 
in the coupled channels of tt~p, Tr°n, r\n, A' + £~, K°Y,°, 
K°A for the neutron resonance (with neutral charge) and 
tt°p, tt + u, rjp, K + I] a , K°H + , K + A for the proton reso- 
nance (with +1 charge). The scattering amp litude for the 
iV(1535) resonance is described in Ref.[17| by means of 
the Bethe-Salpeter equation for meson baryon scattering 
given by 



T = V + VGT 



(18) 



Based on the N/D method and the dispersion relation 
Q, this integral scattering equation can be reduced to a 
simple algebraic equation 



T=(1-VG)~ 1 V 



(19) 



where the matrix V is the s-wave meson-baryon inter- 
action given by the lowest order of the chiral perturba- 
tion theory, which is the Weinberg- Tomozawa interac- 
tion, given by 



u, 



M t + E Mi+E' 



2Mi 



2M, 



(20) 



with the channel indices i,j, the baryon mass M, the me- 
son mass m, the meson decay constant / and the center 
of mass energy y/s. The coefficient CV, is the coupling 
strength of the meson and baryon, which is determined 
by the SU(3) group structure of the channel. The diago- 
nal matrix G is the meson baryon loop function given in 
terms of the meson and baryon propagators by 



M 



1 



(2tt) 4 E(q) ^ - E(q) 
1 

v 

(P — q) 2 — m 2 + it 



(21) 



with the total energy P = (-^,0,0,0) in the center of 
mass frame. For the baryon propagator we use the non- 
relativistic form and neglect the negative energy prop- 
agation. The loop function should be regularized with 
proper schemes. In the practical calculation, we take di- 
mensional rcgularization by using a covariant form of the 
positive energy part of the baryon propagator, 



M Y> r u r {q)u r (q) 2M T, r u r (q )u r (q ) 



E{q) q° - E{q) + it 



M 2 



ze 



(22) 
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In dimensional rcgularization, the loop function in each 
channel i is given by the following analytic expression: 



G, 



d 4 q 



2Mi 



1 



(2tt) 4 q 2 



M 



ie (P — q) 2 — m\ 



167r2 } a^)+\n-t 



2s 



+ ^= [ln(s - {M 2 - m 2 ) + 2q t ^) 
V s 

+ ln(s + (M 2 - mf) + 2q lV Cs) 
- ln(-s + (M 2 - mf) + 2q^s) 
-ln(-s-(M 2 -m 2 ) + 2q l ^- S ) 



(23) 



where q~i is the 3-momentum of the meson or baryon in 
the center of mass frame, p is the scale of dimensional 
rcgularization and cii(p) are subtraction constants, which 
are determined by a fit to the 6*11 and 531 partial waves 
of ttN scattering [llj- Once these constants are fixed to 
the ttN scattering data, the amplitudes involving pho- 
tons can be predicted without introducing any new free 
parameters. 

It should be emphasized that the subtraction constants 
cii ( u) are different for different channels i in the model 
of [17| ■ This is unlike the case of KN scattering and the 
A(1405) resonance, where all the subtraction constants 
in the different channels are approximately equal and of 
natural size according to Ref. |5[ . The need for different 
subtraction constants in the case of ttN scattering and 
the -/V*(1535) resonance has been interpreted recently 
[39I ] as a clear indication that the A r *(1535) contains a 
mixture of a genuine quark state apart from the meson 
baryon components. This conclusion has been reached by 
following an alternative method in which the subtraction 
constants have been chosen of natural order, and approx- 
imately equal, and a CDD pole is included which would 
give us an indication that extra components to the meson 
baryon ones are needed in the iV*(1535) wave function. 
The study of Ref. [39[ clearly indicates that the effect of 
the CDD pole is negligible for the A(1405) resonance but 
relevant for the the case of the TV* (1535). We shall see 
that our approach, based on the meson baryon compo- 
nents exclusively, provides a fair description of data, but 
some remaining discrepancies indirectly hint to the need 
of extra components in the wave function. 

The resulting amplitudes T y from Eq. (fT9|) can be an- 
alytically continued to the complex plane of the scatter- 
ing energy s 1 / 2 . The amplitude has a pole in the complex 
plane that is identified with the resonance, and the cou- 
pling strengths gi of the resonance to the meson-baryon 
channels is determined by the residues of the pole: 



9i9j 



TV 

BG 



(24) 



\/s - M R + iT R /2 
where y/s is the cm. energy of the meson-baryon system 



TABLE I: Complex coupling constants gi of n* to the meson- 
baryon channels. 



7T p 


7r°n 


r\n 


0.557 + 0.325i 


-0.387 - 0.238i 


-1.45 + 0.435i 


K+T,- 


K°X° 


K°A 


2.20 - 0.171i 


-1.56 + 0.115i 


1.39 - 0.0825i 



TABLE II: Coupling constants gt of p* to the meson-baryon 
channels. 



7T°p 


7r + n 


VP 


0.397 + 0.222i 


0.555 + 0.322; 


-1.47 + 0.432* 


K+T, 


K°Z+ 


K+A 


1.56 - 0.133i 


2.21 - 0.183i 


1.37 - 0.100i 



and TgQ is an amplitude for the nonresonant contribu- 
tions. The pole positions of the resonance are obtained 

as 



yfs = 1537 - 37i [ MeV] 
for the n* (neutral charge) and 

Vs = 1532 - 37i [ MeV] 



(25) 



(26) 



for the p* (+1 charge). The values of the coupling con- 
stants are listed in Table |TJ The coupling constants gi 
characterize the structure of the N*. The empirical evi- 
dence of larger coupling of the iV(1535) to r)N than that 
to -kN is reproduced in this model. In addition, the 
couplings to the T,K and AK channels are also large. 
This implies that the AT(1535) has large components of 
strangeness. 

In the meson-baryon picture of the iV(1535) resonance, 
the photoproduction of the resonance from the nucleon 
is formulated through the photon couplings to the me- 
son and baryon components of the A^*(1535). Photon 
couplings and gauge invariance in the case of chiral uni- 
tary amplitudes are discussed in Ref. [ill l4ll. |42| . Here, 
we follow an approach similar to the one developed in 
Refs.[23l. l2~il, [25jj for real photons, extending it to virtual 
ones. Feynman diagrams to the transition form factors 
at one-loop level are shown in Fig. [2] In the loops, all 
possible octet mesons and baryons contribute, namely, 
tt~P, ir°n, r/n, K + Y,~ , K°YP, K°A for the neutron reso- 
nance (neutral charge) and 7r°p, 7r + n, r/p, A" + I] , AT°£ + , 
AT + A for the proton resonance (+1 charge). We sum up 
all the contributions to the transition amplitudes. In dia- 
grams (a) and (b), the photon attaches to the meson and 
baryon in the loop, respectively. Diagram (c) has the 
Kroll-Ruderman coupling which is the contact interac- 
tion of the photon, meson and baryon. Diagrams (d) and 
(e) have to be taken into account to keep gauge invari- 
ance. It seems that these diagrams contain a forbidden 
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P 



(a) 



(d) 




(e) 



FIG. 2: Feynman diagrams for the transition form factor of 
iV(1535) at one loop level. The solid, dashed, wavy and dou- 
ble lines denote octet baryons, mesons, photon and iV(1535), 
respectively. 



transition from the l/2 + state of the proton to the l/2~ 
state of the iV(1535), but negative energy propagation of 
the intermediate baryons in these diagrams is possible, 
having the opposite parity to the positive energy prop- 
agation. The positive energy part in motion also mixes 
different parity states through the different partial waves. 

We calculate the transition amplitudes both in a non 
relativistic and relativistic formulations. The momentum 
of the baryon is small enough to describe the transition 
amplitudes in nonrelativistic formulation. In addition, 
as we have already mentioned above, in the construc- 
tion of the iV(1535) in the meson-baryon scattering [It) ]. 
we have used the nonrelativistic formulation in the ele- 
mentary vertex and the baryon propagators as seen in 
Eqs. (|2"D|) and (|2"3"|) . Therefore, to keep consistency of the 
calculation of the photon couplings with the construc- 
tion of the N* resonance, the nonrelativistic calculation 
is preferable. One should mention that the fact, that 
one fits subtraction constants to data, largely washes out 
relativistic effects from the use of Eq. (f2"Tj) or (f2"3"| in the 
G function in a fair range of energies around the fitted 
point. Nevertheless, it is somewhat complicated to prove 
gauge invariance in the nonrelativistic formalism, since 
we need to take into account all the possible diagrams in- 
cluding negative energy contributions, which are referred 
to as Z-diagrams. To avoid this complication, we will 
also perform the calculation of the amplitudes in rela- 
tivistic formulation, in which the negative energy contri- 
butions are automatically counted without introducing 
the Z-diagrams, and we shall show that the relativistic 
calculation is exactly gauge invariant. This guarantees 
that each term in the 1/M expansion is gauge invariant. 
Exploiting this fact, in the nonrelativistic framework we 
calculate diagrams for leading amplitudes relying upon 
gauge invariance and show that the next-to-leading terms 
are relatively small. 

The basic interactions of the mesons and baryons are 



given by the following chiral Lagrangian: 
D 



MBB 



V2f 
F 

V2f 



Tt[B W {d^,B}] 
Tr[B^ 5 [d^,B}] (27) 



with the meson and baryon fields, <I> and B, defined by 



D 



( -4=7T° + ^-77 

£- 



7T 
^-0 



V 



V2" V& 

K° 
£+ 

"0 



K+ \ 
K° 



(28) 



(29) 



The MBB couplings from these Lagrangian are given 
by g\/{2f) with the corresponding axial vector coupling 
g\ and the meson decay constant /. The axial vector 
couplings arc given in terms of the D and F parameters in 
the Lagrangian (|2"7| as listed in Table HTT1 For the meson 
decay constant, we use / = 93 McV for all channels in 
our calculation. The values of D and F for the axial 
vector couplings are taken from Ref-JU 



D = 0.85 ±0.06 



as 



F = 0.52 ±0.04 



(30) 



These values are determined by the experimental data 
of the hyperon axial vector couplings, neglecting higher- 
order corrections. The photon couplings to mesons and 
baryons are given by the gauge couplings: 

C lB = -eTr [B^[Q ch , B]] A» (31) 
C jM =ieTr[d fi <S>[Q ch ,$]]A» (32) 

with the charge matrix Q c h = diag(|, — |, — -|). The 
Kroll-Ruderman terms, the jMBB couplings, are ob- 
tained by replacing the derivative acting on the me- 
son fields, 9 p 4>, with the covariant derivative D^<& = 
+ ieA^lQch, $] in the interaction Lagrangian (|2T|) to 
realize the gauge invariance. The Kroll-Ruderman terms 
are proportional to the meson charge Qm- For the cou- 
plings of the meson-baryon to the N* resonance, we take 
a Lorentz scalar form representing the s-wave nature and 
the coupling strengths are taken from the chiral unitary 
approach as given in Tables |U [Til 



B. Nonrelativistic formulation 

In the nonrelativistic formulation, the leading terms of 
the 1/M expansions are the diagrams (a) and (c) in FigFJ] 
Diagram (b) is found to be of next-to-leading order due to 
the 1/M factor in the ^BB coupling. In the CM frame of 
the TV*, which we take for the nonrelativistic calculation, 
diagram (d) vanishes, since there is a direct transition 
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TABLE III: The axial vector coupling g\ for each channel. 
The values of D and F are given in Eq. (|30[) . 



channel 


np-K 


nmr° 


nnrj 


9a 


V2(D + F) 


-D-F 


^(-D + 3F) 






nX°K° 


nAK° 




V2{D - F) 


-D + F 


~^(D + 3F) 


channel 


ppiT 


PWK~ 


PPV 


9a 


D + F 


y/2(D + F) 


75 (-D + 3F) 




pE°K+ 


pE+K° 


pAK+ 




D-F 


V2(D - F) 


-js(D + 3F) 



\ N* 



i2Mu(p)u(p) 
p 2 - M 2 + le 



/ 9a 



~\ 2/ 



k ■ a 



V 



i-Qhiip + p') ■ f- 



given in Eq. (fl"4)) . Since the helicity amplitudes can be 
expressed by the A^2 R and M^ K in Eqs. {HJ) and $16]), 
we will calculate only these two amplitudes. The ampli- 
tudes A^^ R and A^^ R remain finite even with one loop 
integration. In fact, .M^ R does have divergence in the 
loop calculation, which should cancel with divergences 
coming from the other diagrams thanks to gauge invari- 
ancc. Here we do not confirm the cancellation of the di- 
vergences, since we later show the complete cancellation 
in the relativistic formulation. 

The Fcynman rules for the nonrclativistic couplings are 
summarized in Fig. [3] In the figure, e M denotes the pho- 
ton polarization, and a is a Lorentz covariant form of the 
spin matrix, — (0, a). g l A stands for the axial vector 
coupling constants of the baryons to the corresponding 
meson. The values of g\ for each channel are given in 
Table IIIII g l N , is the coupling strength of the N* to the 
meson-baryon channel i. The values of g l N » are listed in 
Table U and |TTJ For the baryon propagator we use the 
covariant form from Eq. ([22]) . The jBB coupling is used 
for the calculation of diagram (b), which is not taken 
into account in our final result of the nonrelativistic cal- 
culation. But in order to confirm that sub-leading terms 
from the l/M expansion are negligibly small, we have 
calculated the diagram (b) in the nonrelativistic calcu- 
lation. The ryBB vertex is obtained by a nonrelativistic 
reduction of the interaction Lagrangian (|3ip as 



FIG. 3: Nonrelativistic Feynman rules for the propagator and 
the elementary vertices. The solid, dashed, wavy and double 
lines denote octet baryons, mesons, photon and iV(1535), re- 
spectively. M and m denote the baryon and meson masses, 
respectively. Qb and Qm are the charges of the baryon and 
meson. 



of l/2 + to 1/2 - . Diagram (e) has some contribution in 
this frame, but it is also found to be of next-to-leading 
order, since the contribution to diagram (e) is confirmed 
to vanish in the large M limit. Indeed, if one neglects the 
kinetic energy term of the baryon propagators in the loop 
(p 2 /2M), the loop function vanishes. It also vanishes in 
the rest frame of the nucleon since it again involves a 
1/2+ to 1/2" transition. 

We will obtain each component of the diagram (a) and 
(c) in the decomposition in terms of the Lorentz structure 



£*• (p+p') 



2M 

£■ (P + P r ) 
2M 



X 



X (33) 
(34) 



where we have used the fact that the baryon kinetic en- 
ergies are small in the nonrelativistic kinematics, p° ~ 
p°' ~ M, in the last expression. In Eqs. (J33J) and l[53]) 
we have neglected the magnetic term that behaves like 
(<7 x k)/2M which has one power less in the loop variable. 
In Sec. IV Bl we shall estimate the contributions from the 
convection current and the magnetic terms. In Eqs. (|33[) 
and (|34"1) , Qb is the baryon charge such that it is e for 
the proton with e 2 /(47r) = a ~ 1/137. 

Let us start with diagram (a). Applying the Feynman 
rules shown in Figj3] the amplitude (—it = J ■ e) for 
channel i is calculated as 



it 



d 4 q 
iQ M A 



d A q 



i2Ma 



(p - q y - Mf 



(q - k) ■ a (2q - k) ■ e 



-iQ M )(2q-k) 



[q-kf 



171; 



(2tt) 4 ((P - q) 2 - Mf + ie){q 2 - mf + ie)((q - kf 



■it) 



(35) 



where the coefficient Ai is defined by We use the Feynman parameterization of the integral 

1 

(a + (b — a)x + (c — b)y) 



Ai= g^ N M i ± = 2 fdx Tdy * s . 

J abc In . r 
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Then, using the integral variable q' , such that q = q' + linear terms of q in the denominator and obtain 
P(l — x) + ky and renaming q' as g, we eliminate the 



if a = iQuAil I dx / dy 
Jo Jo 



d 4 q {q+{y- l)fc) • cr (2g + (2y - l)k + 2(1 - x)P) ■ e 



(2tt)4 



(38) 



r 



where we use P ■ a = in the CM frame and S 1 * is defined 
by 

S l a = 2P- k(l - x)y - P 2 x{l -x)- k 2 y(l - y) 

+M?(l-x)+m?x . (39) 



In Eq. (|38p . even powers of q give contributions after per- 
forming the integration. The q^q v term in the numera- 
tor which contributes to the A4^ R is divergent, while the 
terms with Oth power of q remain finite and contribute 
to the A^^ R and A^^ R amplitudes. Finally after the in- 
tegration, we obtain the an( i .Mlp 1 components for 
the channel i as 

i(NR) _ QmM 



M 



2u 



(Air) 2 

where we have used 

A 1 



dx 



D„ It 



dx I dy 



(y-l)(2y-l) 
Si - ie 



(41) 



1 / 1 



(2tt) 4 (q 2 - Sf {±it) 2 2\S 



(42) 



In a similar way we evaluate the contribution from di- 
agram (b) which, as we mentioned, is of order 1/M of 
the previous ones and we obtain 



M 



s(NR) 
2b 



M 



i(NR) _ 
36 — 



QbA 

' (4tt) 2 
' (4tt) 2 



dx 



dx 



dy 



2y(l~x) 
Si-ie 



(43) 



IC 



with 



St = 2P- k(l - x)y - P 2 x(l - x) - k 2 y(l - y) 

+m 2 {l-x) + M 2 x . (44) 



For the diagram (c), the amplitude has only the .Mf R 
component as seen in 



= iQ M A 



d 4 q 



i2M 



(2^)4 (P - g) 2 
d 4 q 



M 2 + ie 
e • a 



(2tt) 4 ((P - g) 2 - M 2 + ie){q 2 - m 2 + ie) 



(45) 



Hence we do not need to perform further calculation for 
this amplitude. 

Finally the helicity amplitudes in the nonrelativistic 
formulation arc obtained by summing up all the chan- 
nels and substituting the amplitudes (|4"0"]) and (T4"Tj) in 
Eqs. (fT5 | and (fTe)) . 



C. Relativistic formulation 

In this subsection, we calculate the transition ampli- 
tudes in relativistic formulation at one loop level. One of 
our purposes for the relativistic calculation is to confirm 



gauge invariance of our formulation. Without the 1/M 
expansion, which has been performed in the nonrelativis- 
tic calculation, all the diagrams shown in Fig|2] should be 
calculated to make the amplitudes gauge-invariant at one 
loop level. Each diagram has divergence from the loop 
integral. It will be found that the divergence appears 
only in the Mi term. After testing gauge invariance by 
summing up all the diagrams, the amplitudes should be 
finite without any rcgularization, since gauge invariance 
assures cancellation of the divergences coming from each 
diagram. In order to check this cancellation of the diver- 
gences, we calculate the M.\ terms from all the diagrams. 
To isolate the divergent parts of the amplitudes, we ex- 
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i($ + M) i 
p p 2 — M 2 + ie p p 2 — m 2 + it 




■ e 



ploit dimensional rcgularization, which respects gauge- 
invariance. The calculations are done in d-dimension, 
and then we expand d around d = 4 in terms of e given 
by d = 4 — 2e. We also calculate the finite and 
to obtain the helicity amplitudes in the relativistic for- 
mulation. We will find again that only diagrams (a) and 
(b) contribute to the ftA\ and M.2 amplitudes. 



Let us start with diagram (a). Using the Feynman 
FIG. 4: Same as Fig[3]for the relativistic formulation. Q E TU \ CS shown in Figd the amplitude of the diagram (a) 
stands for the external N* charge. for the channcl l is givcn by 



- iTl = 



f d d q 

IQmB; 



-igi)i 
d d q 



f-i+ M. t 



{P-qf 



-4 + Mi 



le 



9a 
2/ 



7n «-fO-nsi 



k) 



1 



{q-kf 



it q z — m- + le 



{2ttY {{P - q) 2 - M 2 )((q - k) 2 - m 2 )(q 2 - m?) 



-«Qm)(2(? - k) ■ e 
(46) 



where we define Bi = {g l A g l N , /2f). We can write ex- are found as 
pressions for all the other diagrams and then an explicit 
calculation shows that by substituting e M by k^ one ob- 
tains an exact cancellation of the terms, hence passing 
the ordinary test of gauge invariance. The coupling of 
the photon to all lines and vertices in the meson baryon 
loop diagrams guarantees gauge invariance, as is com- 
monly known [HI, [H[ . j 

We now come back to the amplitude T a of Eq. ([46]) and, 
using the Feynman parameter integral (|37p . shifting the 
integral variable q to q = q' + P(l — x) + ky and renaming 
q' by q, we find the integrand written as a function of q 2 

as 



iT* = iQ M Bi2 



d;r 



dy 



\q 2 C a + D a 



(2n) d (g2 _ ft + le f 



(47) 



D 



a 3 



fit 

c\ 



a2 



c 



a3 



(M N + Mi)(l + %) 
2(3x-2) + xe 
2(1 - 3y) - ye 



D* al =0 



2(1 - x) [{x - y)yk 2 - x{x - y)M%. 
+y{l-x)M 2 N + {y-l)MiM N 
-(x-y)M N *(M N + Mi)} 
{2y - 1) [(a - y)yk 2 - x(x - y)M 2 N , 
+y(l - x)M 2 N + (y- hU.W.v 
-{x-y)M N ,{M N + M i )} 



(50) 
(51) 
(52) 



(53) 



(54) 



(55) 



where S a is defined in Eq. (|39|) and the coefficients C a 
and D a are defined in terms of the Lorentz components 
by 

C a = C al ?j 5 + C a2 P -675 + C a3 fc-e 75 (48) 
D a = + D a2 P ■ e 75 + D a3 k ■ £75 . (49) 

After some algebra, noting that g£ = d = 4 — 2s in the 
d-dimcnsional calculation, the coefficients C a i and D a i 



After the q integration, the divergent term is found only 
in the M.\ amplitude as 



QuBj 1 M 
' (4tt) 2 e ' 



N 



Mi 



ii5 



(56) 



The divergent terms in the M.i and M. 3 vanish due to 
/„ dx /* dy(3x- - 2) = and da; /* <%(1 - 3y) = 0. To 
calculate the divergent term we have used the following 
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formula: 



/ 



d d q 



4-2er(e) (1 



The finite parts in the M.2 and M3 terms are obtained 
as 



{2n) d (q 2 - S*) 3 (47r) 2 - £ 2 r(3) \S / 
j^y Q - log5 - \ ~ 7 + log 4tt + 0(e)) (57) 
I 



2« 



(47T)" J0 



2 

^Sa" (47r)2 



/ c?y 
'0 



2(3x- 2)logS* -x 



D„ 



dx I dy 
Jo 



2(1- 3y) log + 



(58) 
(59) 



where the first two terms in the integrands are from the finite parts of the divergent integrals and the last terms come 
from the finite integrals. 

Next, let us move to the calculation of the diagram (b). In a similar way. we obtain the finite M.2 and M3 
amplitudes as 



with S b defined in Eq. 



D 



M\ b = 
M\ b = 

and 

1,2 



"(47)2" 

QbB, 
(4tt) 2 



dx / dy 



D 



2(3y-l)log^-y + l + ^ 

J, 



d.v 



dy 



-2(3» -1) log Sj + y-1 



b 

D b3 



2(y(l - y)(x - y)k 2 - y(l - x)(x - y)M 2 N , - y 2 (l ~ x)M% 
+xyM N Mi + (x- y)((l - x)M N , + .\/,!i.\/ v + M*)) 
D l b3 = 2(-y(l - y)(x - y)k 2 + y(l - x)(x - y)M 2 N , + y 2 (l - x)M\ 
-xyM N Mi + y(x + y)(M N + Mi){M N * + M N )) 

I 



(60) 
(61) 



(62) 
(63) 



We have the divergent term in the M.\ term as 
Q B Bi 1 Mat - M N , - M, 



(4tt) 2 e 



^75 



(64) 



The amplitudes for diagrams (c), (d) and (e) have only 
the A4\ components, which we do not use for the transi- 
tion amplitudes. Only the divergent terms are necessary 
for the present arguments. The divergent terms are found 



according to 



A=a 



(4tt) 2 







{Qe — Qb — Qm) 



M N * - M N + Mi 1 



£ 

(68) 



due to the charge conservation Qe = Qb + Qm- The 
cancellation takes place in each channel of the loop. 



Co 



Q M Bi 1 Mat. + 2M, 
-jA^fS 2 ^ 

Q E B t 1 (M%. - 2M 2 

(4tt) 2 £ 
Q E B l 1 (M 
(4tt) 2 e 



2m 2 + M N *Mi) 



2(M N + M N . 
2M 2 - 2m 2 - M N M 



2{M N , + M N ) 



(65) 
(66) 

-Hs ■ (67) 



At the end, collecting all the divergent terms of dia- 
grams (a) to (e), we find that the divergent terms cancel 



IV. RESULTS 

In this section we show our results for the helicity am- 
plitudes, A1/2 and S1/2, of the iV(1535) dynamically gen- 
erated in mcson-baryon scattering. In Fig|5l we show our 
result for the A 1 / 2 amplitude of the proton resonance cal- 
culated in the nonrelativistic formulation [cf. Eqs. (fl~5l 
HOI SU)] with the CM energy W = 1535 McV. In the 
present calculation, we multiply the amplitudes obtained 
in the former section by the electromagnetic form factors 
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dependence comes only from the loop calculation per- 
formed in the previous section. In this case, the helicity 
amplitude increases as Q 2 increases. The inclusion of the 
form factors introduces a decreasing function of Q 2 which 
leads to a Q 2 dependence of the helicity amplitude in fair 
agreement with the experimental observation, although 
it falls faster than experiment since at Q 2 = we need a 
renormalization factor of 1.45 to reach the data while at 
Q 2 — 1 GeV 2 we need a factor of 2.15. 

The absolute magnitude of our helicity amplitude looks 
underestimated if one compares our result directly with 
the experimental data shown in the figure. But it should 
be noted that extraction of the helicity amplitude from 
the experimental observables of the jp — > rm reaction is 
performed by using the following formula [IH, EH, HH, [HH : 



FIG. 5: Modulus of the A\/2 helicity amplitude for the pro- 
ton resonance as a function of Q 2 with W = 1535 MeV cal- 
culated in the nonrelativistic formulation. The solid (dotted) 
line shows the A^ 2 amplitude with (without) the form fac- 
tor of the meson inside the loops given in Eq. (|69p . Marks 
with error bars are experimental data normalized by the N* 
full width rV = 150 MeV and the iV* -> r/N branching 
ratio b n = 0.55. Filled triangles and circles are results of 
the CLAS collaboration taken from Refs. 133] an d |45| , r espec- 
tively. Open circles show results of Refs. [4(1 l47l . [48, 49]. The 
values are taken from Ref.Bol. 



of the mesons or baryons to which the photon couples. 
The form factors of the meson and baryons components 
of the resonance, together with the intrinsic Q 2 structure 
of the loops are responsible here for the Q 2 dependence 
of the helicity transition form factors. For the mesons 
and baryons form factors we take monolopc form factors 
consistent with the values for the radii of the mesons. We 
take 



F(Q 2 ) 



A 2 



A 2 + Q 2 



with 



K = 0.727 [GeV] 
A K = 0.828 [GeV] 



(69) 



(70) 
(71) 



which correspond to (r 2 ) = 0.44 fm 2 and (r 2 ) = 0.34 fm 2 
for the pion and the kaon respectively [H(3, EI H3] • For 
the baryon form factor, we take the same form as for the 
corresponding meson to keep gauge invariance. 

In FigJSJ we show our result for the A^, 2 amplitude 
of the proton resonance calculated in the nonrelativis- 
tic formulation together with various experimental data. 
The CM energy is taken as W = 1535 MeV. Let us 
first discuss the Q 2 dependence of the helicity ampli- 
tude. The solid line denotes the modulus of the cal- 
culated amplitude multiplied by the meson form factor 
given in Eq. ([B"9"|) , while the dotted line shows the results 
without the meson form factor, which means that the Q 2 



A 1/2 (Q 2 ) 



IWT 



N" 



2m p b ri 



a(W,Q 2 ) 



(72) 



with a N* full width T N , , a N* — > -qN branching ratio 
b v , a resonance part of the total cross section a(W,Q 2 ), 
the CM energy W and the proton mass m p . To obtain 
this relation, one assumes that the cross section is domi- 
nated by the single 7V(1535) resonance and that the Si/% 
amplitude is small. For the experimental data shown in 
Figj5[ the amplitudes are normalized by T^r. = 150 MeV 
and b v — 0.55 (45|, which are obtained in a global fit of 
the cross section with the Breit-Wigner amplitude. On 
the other hand, the TV* width obtained in the present ap- 
proach is Tat. ~ 74 MeV for the p* as obtained from the 
pole position shown in Eq. ([26[) . in which the half width 
is given by the imaginary part. The branching ratio b v 
in this approach has been reported as b v ~ 70% fl7j |. 
This normalization difference would give us a factor 1.6 
reduction in the normalization with respect to the data 
shown in Fig[5j Similarly, should one use in the experi- 
mental analysis a 7V(1535) width of the order of 90 MeV 



as found at BES [57 
MAID2007 analysis 



or the 100 MeV quoted in the last 
58| , the results obtained would be in 
much better agreement with the theoretical results. We 
will come back to the discussion on the normalization 
later on when discussing the photoproduction cross sec- 
tion in the present approach. Let us note that the value 
obtained here at Q 2 = of A p 1/2 = 64.88-10~ 3 GeV" 1 / 2 is 
in excellent agreement with the most recent MAID 2007 
analysis reported in Ref.g^ of 66 • 10" 3 GcV _1/2 . 

The Sj^ 2 amplitude calculated in the nonrelativistic 
formulation is plotted in Fig[6] together with experimen- 
tal data. Although the modulus of Si/ 2 is plotted in the 
figure, the ratio of S'1/2 to A1/2 is nearly real and negative 
in agreement with experiment, where an implicit phase 
convention is taken that renders A^, 2 real and positive. 
Here we also show the effect of the meson form factor. 

As mentioned before, diagram (b), in which the pho- 
ton couples to the baryon in the loop, gives sub-leading 
contributions in the l/M expansion. This can be seen 
in Figs [HI M an d [TTJ] In these figures, the amplitudes 
with diagram (a) only (solid line) arc almost equivalent 
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to those with both diagrams (a) and (b) (dashed line), 
and the contributions from diagram (b) (dotted line) are 
smaller, in the present case, than typical corrections of 
20-30 percent for the 1/M terms. Therefore, the helic- 
ity amplitudes in the nonrelativistic formulation around 
these energies are basically given by diagram (a). We also 
plot the pion and kaon contributions separately. The fig- 
ure shows that the pion contribution (dot-dashed line) is 
comparable with the kaon contribution (two-dotted line). 
This implies that the strange component is important in 
the structure of the iV(1535). 

Next, we report on results for the hclicity amplitudes 
of the neutron. The n/p ratios of the hclicity amplitudes, 
A™/ 2 /A^, 2 and Sy 2 /S P ^ 2 , are plotted in Fig[Tl]as a func- 
tion of Q 2 . For a real photon at Q 2 = we obtain the 
ratio —0.79 + O.lli, which is almost a real value, and its 
modulus 0.80. A multipole analysis [59| using the inclu- 
sive experimental data of Ref. [46[ gives the negative sign 
value -0.84±0.15 for A™ /2 /A p 1/2 . Values of \A™ /2 \/\A p 1/2 \ 
which are extracted from the ratio of the eta photopro- 
duction cross sections, a n /a v , are reported as 0.82 ±0.04 
in Ref.[g3| and 0.819 ± 0.018 in Ref.0. The result ob- 
tained in our approach agrees with the experimental data 
in both sign and magnitude. This comparison is free from 
the normalization uncertainty of Eq. (|72[) . 

The values of the A X j 2 hclicity amplitude at Q 2 = are 
summarized in Table IIVI The phases of the amplitudes 
are set so that the A 1 / 2 hclicity amplitude for p* has 
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FIG. 7: (Color Online) A^ 2 helicity amplitudes for the pro- 
ton calculated in the nonrelativistic formulation as a function 
of Q 2 . The upper, middle and lower panels are respectively 
the modulus, real parts and imaginary parts of the ampli- 
tudes. The phases of the amplitudes are set so that the -A^ 2 

amplitude has a real and positive value at Q 2 = 0. The 
solid lines show the calculation with diagram (a) (meson pole 
term). The dashed and dotted lines stand for the calculations 
of sum of diagrams (a) and (b) and diagram (b) only, respec- 
tively. The dot-dashed and two-dotted lines denote the pion 
and kaon contribution, respectively. 



FIG. 6: Modulus of the 5*1/2 helicity amplitude for the proton 
resonance as a function of Q 2 with W = 1535 MeV. The solid 
(dotted) line shows the S^ 2 amplitude with (without) the 
form factor of the meson inside the loops given in Eq. (J69}. 
The sign of this amplitude relative to A^ 2 is negative, both 
in experiment and theory. Solid triangles up (down) show the 
results from a combined analysis of n (r/) electroproduction 
data J27|, HH . The solid squares are from single-Q 2 fits from 
Ref. [H3]. The empty triangles up are taken from Ref. 
The other data (empty circle and empty triangle down) are 
from Ref. [H. 



a real and positive value at Q 2 = 0. The ratios of the 
helicity amplitudes to that of p* are also shown in the 
table. We also show the hclicity amplitudes in the isospin 
decomposition: 

4% = ^(4/2+^i/ 2 ) (73) 
^2 = l(Al /2 -A? /2 ). (74) 

In the nonrelativistic calculation we find that the value 
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FIG. 8: (Color Online) S^ 2 helicity amplitudes for the proton 
calculated in the nonrelativistic formulation as a function of 
Q 2 . Same as Fig[7] 



of the isoscalar component is much smaller than that 
of the isovector, which is consistent with experimental 
observation. The ratio of the isoscalar component to the 
p* amplitude is \Af /2 /A p 1/2 \ = 0.12 in our calculation. 

while in experiments it is found to be A l ^ /2 /A p 1/2 = 0.09± 
0.02 in Ref.fHi and 0.09 ± 0.01 in Ref.[6l[. 

It is also interesting to compare the values of our p* 
and n* helicity amplitudes A 1 / 2 at Q 2 = with those 
of the PDG [62]. We obtain 0.065 GeV" 1 / 2 and -0.052 
GeV -1 / 2 for the p* and n*, respectively, versus the val- 
ues quoted in the PDG, which include uncertainties from 
the compilation of data of several analyses, 0.090 ± 0.030 
GeV" 1 / 2 for the p* and -0.046 ± 0.027 GeV" 1 / 2 for the 
n* . As one can see, the agreement, within uncertainties, 
is good. 



FIG. 9: (Color Online) A™/ 2 helicity amplitudes for the neu- 
tron calculated in the nonrelativistic formulation as a function 
of Q 2 . The upper, middle and lower panels are respectively 
the modules, real parts and imaginary parts of the amplitudes. 
The phases of the amplitudes are set so that the A^ 2 ampli- 
tude has a real and positive value at Q 2 = 0. The solid shows 
the calculation with the diagram (a) (meson pole term) . The 
dashed and dotted lines stands for the calculations of sum of 
diagrams (a) and (b) and diagram (b) only, respectively. The 
dot-dashed and two-dotted lines denote the pion and kaon 
contribution, respectively. 



V. DISCUSSION 

A. Photoproduction of the eta meson 

In this section, we investigate the photoproduction of 
the eta meson close to threshold energies in order to dis- 
cuss the normalization of the helicity amplitude in the 
present approach. In the calculation of the helicity am- 
plitude in Sec. IIII Bl we have separated out the 7V(1535) 
resonance contributions from the scattering amplitudes 
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FIG. 10: (Color Online) Si/ 2 helicity amplitudes for the neu- 
tron calculated in the nonrelativistic formulation as a function 
of Q 2 . Same as FigH 



in which the 7V(1535) is dynamically generated, by set- 
ting the CM energy as the resonance energy and multi- 
plying the coupling strengths of the 7V(1535) resonance 
to each channel, g l N » , by the loop functions. For our pur- 
pose of calculating the r\ photoproduction cross section, 
we replace the coupling strengths, g l N , , by the MB — > rjp 

scattering amplitudes, t^p, obtained by the chiral unitary 
approach [T3|, where i denotes the initial meson baryon 
channel. The evaluation of these amplitudes is sketched 
here in section Mil and done in detail in Ref. [I?) ■ 

Following the above prescription for the eta photopro- 
duction amplitudes, we obtain the cross section of the 
photoproduction as 



M 2 ^n_\ + i2 
Airs ky 1 7P ^" pl 



(75) 
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FIG. 11: The n/p ratios of the helicity amplitudes with W = 
1535 MeV. The upper and lower panels show the np ratios of 
the A 1 / 2 and Si/ 2 amplitudes, respectively. The open square 
shows the value given in Ref. [591] . 



Nonrelativistic calculation 




A1/3 


\A 1/2 \ A 1/2 /A p 1/2 


\Ai/2/Al /2 \ 


P* 


64.88 


64.88 




n* 


-51.54 + 7.2H 


52.04 -0.79 + 0.11i 


0.80 


IV 


58.21 - 3.6H 


58.32 0.90 - 0.056* 


0.90 


IS 


6.67 + 3.6H 


7.59 0.10 + 0.056i 


0.12 


Relativistic calculation 




A 1/2 


l^i/2 1 A 1/2 /A p 1/2 


\Ay,2/Al /2 \ 


P* 


46.31 


46.31 




1 1 * 


-55.24 + 28.95i 


62.36 -1.19 + 0.63? 


1.35 


IV 


50.78 - UA7i 


52.80 1.10 - 0.3H 


1.14 


IS 


-4.46 + UA7i 


15.15 -0.10-0.31i 


0.33 



TABLE IV: Values of the A 1/2 helicity amplitudes at Q = 
in units of 10 _3 GeV -1 ^ 2 in the nonrelativistic (upper panel) 
and relativistic (lower panel) calculations. The ratios to the 
A^, 2 are also shown. The phases of the amplitudes are set 
so that the A^, 2 amplitude has a real and positive value at 
Q 2 = 0. IV and IS stand for isovector and isoscalar. 



CM frame and the T-matrix is given by 



o 7p— >r)p 1 



|2 _ 8ml + 8EE' 
\Qrril 



6 .(*) 



(76) 



where fc 7 {k v ) is the photon (rf) three-momentum in the where E (£") are the energies of the incoming (outgoing) 
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FIG. 12: Cross section for photoproduction of the r\ using the 
nonrelativistic formalism for the photon loop. Dots: Analysis 
from Arndt [H (SAID data base). 



proton in the CM frame. To obtain the «M]/ N am- 
plitude from the At^ 1 * and K ' amplitudes calcu- 
lated in the previous section, we use the gauge invariance 
condition given in Eq. (|17p . Actually, the eta photopro- 
duction with a real photon has no contribution from the 
A^3 amplitude, hence only one amplitude needs to 
be evaluated which we choose to be M.^ that shows its 
finitcness immediately. 

In Fig[T21 the total cross section of the present calcu- 
lation (solid line) is plotted together with the data from 
Arndt [63| (dots) . Our result of the eta photoproduction 
cross section provides the right strength around the peak 
of the A r *(1535) resonance but the width of the peak is 
narrower than the experiment as a result of the narrower 
widths of the 7V(1535) resonance obtained by the present 
model. 

In Fig[T3] we show the ratio of the cross sections of rj 
photoproduction on the neutron over that on the pro- 
ton, (7 n /a p , obtained in the nonrelativistic formulation, 
comparing our calculation with experimental data. The 

785 MeV which corresponds 
1535 MeV is found to be 0.53, which is quite 
consistent with the experimental data. 

Although we are only concerned with the vicinity of 
the iV*(1535) resonance, one cannot overlook the appar- 
ent discrepancy of the theory and experiment at photon 
energies above 800 MeV as shown in Fig. [T31 which is 
also shared by the model of Ref. [2| . Only very recently 
have we obtained experimental information that brings 
a new perspective to these discrepancies. Indeed, in Ref. 
[65j the theoretical results of Fig. [TJlare taken and folded 
with the Fermi motion of the nucleons in the deuteron to 
allow a realistic comparison with the experimental data. 
In 65j it is shown that the steep rise of the theoretical 
curve is softened to a curve in between the one of Fig. 
[T3l and a horizontal line starting from 800 MeV. On the 



value of the ratio at E}^ h 
to Eqm 



FIG. 13: Ratio of cross sections of photoproduction on the 
neutron over that on the proton, a n /a p , as a function of the 
photon energy E-, in the laboratory frame. The data are 
taken from Ref. [6l| for the deuteron target and Ref. [U for 
the Helium target. The dashed line is a theoretical calculation 
by Kaiser et al. in Ref. @]. 



other hand, recent results from [66| for jn — > r\n show a 
steady rise starting from E 1 = 900 MeV. These two facts 
together would render the apparent discrepancies into a 
rough qualitative agreement. Let us mention in this re- 
spect that the inclusion of the 7T7rAf channel, although 
only qualitatively considered as shown in the next sec- 
tion, also works in the direction of softening the steep 
rise of Fig. 1131 In any case, we must admit larger the- 
oretical uncertainties at higher energies than around the 
resonance region, also including extra terms considered 
in Ref. [i^] that would become relevant as one moves 
away from the resonance pole. 



B. Higher order couplings 

In the nonrelativistic treatment of the photon loops 
from Fig. [2] in Sec. IIII B[ the magnetic couplings of the 
photon to the baryons have been neglected as they are of 
higher order in the external photon momentum accord- 
ing to k/M. The convection part of the ^BB coupling 
shows a similar p/M suppression, where p is a typical 
loop momentum. These higher order terms have been 
neglected for the sake of consistency with the hadronic 
part of the model: as discussed after Eq. (f2Tj) . only the 
positive energy part of the baryon propagator is taken in 
the evaluation of the MB — ^ MB scattering amplitude. 
In this section, the effects of the higher order jBB cou- 
pling are studied, which can give an idea of theoretical 
uncertainties from these terms. 
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The photon-baryon coupling is given by 



£ 7 bb = -* Qb4 



2M 



a^d^ * (77) 



with a^ llJ = | [7^,7"], the baryon charge Q B and the 
anomalous magnetic moment n given in units of the nu- 
clear magneton yu^r = e/(2Afjv). In the nonrelativis- 
tic reduction of this interaction, only terms up to order 
p/(2M) are considered which leads to the vertex 



it 



iQb F(Q 2 
2M 



e(p + p>) + 



eG M (Q 2 



2M 



N 



e (k x a)/iB 
(78) 



where we have supplied the form factor F(Q 2 ) given in 
Eq. ([55]) and the Sachs form factor G M (Q 2 ) = 1/(1 + 
Q 2 /A 2 M ) 2 , K 2 M = 0.71 GcV 2 and fi B is the baryon mag- 
netic moments in units of /in from the PDG [6^ . We 
use common form factors Gm for the £ and A. 

With the vertex from Eq. ([75]) . diagram (b) from Fig[5] 
can be calculated. For the convection part, the result has 
been already obtained in Eq. (p4"3"]) . The magnetic part of 
diagram (b) is given by 



''mag 



V 



e (a x k) 



(2^) 4 q 2 -M 2 



2M r 



(q-k) 2 -M 2 {P-qf 



(79) 



This expression is finite and not logarithmically divergent 
as the convection part of the ^BB coupling or diagrams 
(a) and (c). The magnetic part is gauge invariant by itself 
as the structure of Eq. (|T9"|) shows. In Eq. ([79"]) . m(M) are 
the masses of the meson (baryon) of channel (i), P 2 = 
s, and g l N , are the coupling strengths to the jV(1535) 
from Tables Q] and [Til For photoproduction, one obtains 
the amplitude T{-yN — > r)N) by replacing g l N , with the 
MB — > Nrj T-matrix as discussed in Sec IV A[ summing 
over all channels i. The axial charges g\ are given in 
Table IIIII We also take into account the magnetic S°A 
transition. We choose a negative ^soa = — 1-61 which 
is the prediction of the quark model [631 while only the 
modulus can be measured [II] . For the photon loops with 
a S°A transition, we use average masses for the baryons. 
For the unknown magnetic moment of the TP, we take 
Ms = 5(^1;+ + My;-- ) = 0.65 which is obtained by the 
SU(3) argument [681 ] and also consistent with the quark 
model. 

Evaluating the loop integral using Fcynman parame- 
ters we obtain 



mag 



with 



ih {i) (k 2 a -e-k-ak-e 



dx I dy 



(») 9n* 9a M i 
16tt 2 / 



It 



(80) 



(81) 



TABLE V: The 1/M contributions from the baryon pole term 
for A 1/2 {Q 2 = 0) [10 _3 GeV _1/2 ]. See text for the different 
cases. 



I 

11 
III 

IV 



A l/2 

64.7 + 4i 
67.7 - 2i 
70.0 - 2.7i 
74.7 



An 

-51.9 + ii 
-42.7- 5. 5i 
-39.9 + 4.6i 
-44.5 + 1.9i 



where is defined in Eq. 



and we used P ■ a = 



in the CM frame. After summing up all channels for the 
photon loop, we obtain the contribution to the helicity 
amplitudes from the magnetic couplings according to 



41/2 



S" 1 / 2 = 

u mag 



2na y/2k 2 ^ , (i) 



'in 



(82) 
(83) 



The jBB coupling of diagram (b) from Fig[2] has a con- 
vection part and a magnetic part (cf. Eq. ([78]) ). and both 
are of order 1/M. Thus, we will treat both parts together 
and compare them to the previous results when only the 
leading order couplings are considered. The latter ap- 
pear in diagrams (a) and (c), while in the nonrelativistic 
calculation diagrams (d) and (e) do not contribute. 

In order to see the effects of the various 1/M terms, 
we have compared them in Tab. [V] The phase is chosen 
in the way that A^ 2 (Q 2 = 0), including all 1/M terms, 
is real and positive. Case (I) is from diagrams (a) and 
(c) of Fig. [2] only, without any jBB couplings. Case (II) 
also includes the convection part of the ^BB coupling 
from diagram (b). Case (III) additionally includes the 
magnetic part. Case (IV) includes, on top of the other 
contributions, the E°A transition magnetic part. 

From the different contributions, we can see that the 
jBB coupling has moderate influence on the results: For 
A^, 2 , the result increases due to the convection part and 
the magnetic part. For A^ 2 , the various 1/M contribu- 
tions lead to a decrease as Tabic IVl shows. As a result, 
the ratio ^1/2/^1/2 = — 0-60 + i0.03 is smaller than the 
value of —0.79 + i 0.11 found in the last section, where 
only the leading couplings were included. 

In Fig. [14] the ratio of 77 photoproduction on the neu- 
tron over that on the proton is shown. The decrease, 
when including the higher order couplings of diagram 
(b), reflects the results from Table IVl that A^ 2 increases 
and A™^ 2 decreases when including the 1/M corrections. 
In the same figure, we show the result when including ad- 
ditional ingredients for the rescattering model from [l7| 
(dotted line). These are the ttttN channel in rescattering 
as well as a form factor for the MB — > MB transitions. 
Coupling the photon to these ingredients is beyond the 
scope of this work; therefore, they are not included in the 
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FIG. 14: (Color online) The ratio a n /a p including the baryon 
pole diagram (b) (dashed line), compared to the result with- 
out diagram (b) (solid line). Also, the result is shown when 
the 7T7rN channel from 11711 is included. The data is from Refs. 
[U,|6l] [see also Ref. [gg]. 

final results. Including these ingredients in the rescat- 
tcring part, the ratio drops for higher photon energies, 
and this gives an idea of theoretical uncertainties coming 
from omitting the ttt:N channel in the present study. For 
other obscrvablcs discussed in this study, the additional 
ingredients from fl7| lead only to very minor changes. 

For consistency with the approach followed for me- 
son baryon scattering, the terms of order 1/M should 
be omitted, and the results of the former section should 
be used to compare with data. The discussion in this sec- 
tion gives us an idea of the uncertainties that one may 
have when considering the 1/M terms. 



C. Result of the relativistic formulation 

Here we briefly discuss our result for the A^ 2 ampli- 
tude of the p* resonance calculated in the relativistic for- 
mulation given in Sec lHI Cl As we have already men- 
tioned, the relativistic formulation is less consistent with 
the model of the 7V(1535) resonance generated dynam- 
ically in the present approach than the nonrelativistic 
formulation. Therefore we rely upon the nonrelativistic 
calculation. However, the Q 2 dependence could be given 
better by the relativistic calculation, particularly if we 
go to values of Q 2 of the order of 1 GeV 2 or above. On 
the other hand, in the relativistic formulation we have 
found the cancellation of divergences coming from each 
diagram. In Fig[15]we show the results for the Ay 2 am- 
plitude in the relativistic calculation shown by the thick 
solid line in comparison with the nonrelativistic calcula- 
tion shown by the dashed line. The relativistic result is a 
bit below the nonrelativistic calculation. We also plot the 
results for the Sf, 2 amplitude obtained in the relativistic 



formulation in Fig. [IBJ The relativistic calculation gives 
a smaller result than the nonrelativistic one, as in the 
case of the , 2 amplitude, but the differences are now 
larger. This reflects the fact that the Sy 2 amplitude is 
more sensitive to small changes of the input and, conse- 
quently, one must accept larger theoretical uncertainties 
in this amplitude. The dispersion of the data seems to 
reflect a similar problem on the experimental side, the 
results proving also rather sensitive to the assumptions 
made in the different analyses. 

Although we prefer the nonrelativistic results of A^ 2 

at Q 2 = 0, as already mentioned, one can take this dif- 
ference as a measure of the theoretical uncertainties. At 
Q 2 = the relativistic result is some 25% below the 
non relativistic one. Taking into account 10% of uncer- 
tainty in the recent MAID2007 analyis [70| of A\ /2 = 

66 ± 7 • 10~ 3 GcV -1 / 2 , one finds good agreement of the 
theory with the MAID2007 results within uncertainties. 

Let us note that in the relativistic calculation the factor 
needed to agree with data at Q 2 = is 1.9 and at Q 2 = 1 
GeV 2 the factor needed is 2.4. It indicates a faster fall 
down than experiment, 25 % lower than exeriment at 
Q 2 = 1 GeV 2 with a curve normalized at Q 2 = 0. This 
compares with 48 % smaller strength than experiment at 
Q 2 = 1 GeV 2 of the nonrelativistic curve, normalized to 
the data at Q 2 =0. This indicates that relativistic effects 
play some role at large Q , as one might think, along 
the line of similar findings in relativistic quark models 

We have also calculated the transition amplitudes us- 
ing N*BM couplings with a derivative of the type j^d^ 
in the relativistic calculation. In this case, we have an 
extra momentum in the N*BM coupling and this mo- 
mentum is included in the loop integral. We also have an 
extra Kroll-Ruderman contact term. As a consequence, 
the cancellation of the divergences is not complete in the 
case of Q 2 > 0, while, for the real photon, that is Q 2 = 0, 
the sum of the amplitudes is finite. From the viewpoint 
of consistency with the model of the N* , one should not 
use the derivative coupling in the N*BM vertex. In the 
unitarization based on the N/D method, we exploit the 
so-called elastic unitarity, in which the interactions V 
are evaluated on the mass shell. In the present case, 
the momenta in the Weinberg- Tomozawa couplings are 
set on the mass shell. Therefore, to maintain consistency 
with this procedure, the momentum in the N*BM should 
have the value on the mass shell and should not be in- 
cluded in the loop integral. This means that a constant 
N*BM coupling is more consistent with the present N* 
model. In any case, just for illustrative purposes, the 
value that we obtain for A^, 2 with the off shell deriva- 
tive coupling is of the order of 90 ■ 10 -3 GeV -1 ' 2 . Our, 
preferred, nonrelativistic result lies in between these two 
illustrative relativistic results. 
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FIG. 15: Modulus of the Ax/2 helicity amplitude for the pro- 
ton resonance as a function of Q 2 with W = 1535 MeV in 
relativistic formulation. The solid (dotted) line shows the 
A^ 2 amplitude calculated in the relativistic (nonrelativistic) 
formulation. The marks are same as in FigO 
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FIG. 16: Modulus of the S1/2 helicity amplitude for the pro- 
ton resonance as a function of Q 2 with W = 1535 MeV in 
relativistic formulation. The solid (dotted) line shows the 
Sy 2 amplitude calculated in the relativistic (nonrelativistic) 
formulation. The marks are same as in Fig[6] 



VI. SUMMARY AND CRITICAL 
OBSERVATIONS 



In this work we have addressed the evaluation of the 
electromagnetic helicity form factors for the electropro- 
duction of the TV* (1535) resonance considered as a dy- 
namically generated resonance. For this purpose the cou- 
pling of the photon to the meson baryon components of 
the N(1535), previously studied within the chiral unitary 
approach to pion nuclcon scattering, was considered. The 
calculations have been done relativistically and nonrcla- 
tivistically, and both of them are found to lead to finite 



results for the transition amplitudes, as well as for eta 
photoproduction which is evaluated simultaneously with 
the same formalism. 

Our study finds interesting results which we summarize 
here: 

The amplitudes were obtained without any free param- 
eters, since the couplings of the resonance to the channels 
have been obtained from a previous study of 7rTV scat- 
tering. The agreement with the A^ 2 amplitude of the 
proton TV* (1535) resonance is fair up to the normaliza- 
tion problem that we have discussed. Indeed, we showed 
that the absolute values of the experimental amplitudes 
were tied to assumptions on the total width of the res- 
onance, which is still far from being a settled issue. We 
also showed that our results for Ai/ 2 at Q 2 = are in 
perfect agreement with the most recent MAID2007 anal- 
ysis of scattering and photoproduction data. The Q 2 
dependence of the transition form factor obtained was 
in fair agreement with the experimental determination, 
although it provided a moderately faster fall down than 
experiment. This result is by no means obvious within 
the picture of a dynamically generated resonance, since 
the Q 2 dependence should be provided by the meson form 
factors and they fall much faster than these experimental 
form factors. Yet, we found that the theory, in the ab- 
sence of the meson form factors, provided a rising func- 
tion of Q 2 , due to the structure of the loops involved, 
which led to a moderate decrease of the TV* (1535) tran- 
sition form factors when the meson form factors were 
considered. 

The results obtained for the S\/ 2 amplitude are also in 
fair agreement with experiment, both in size and the rela- 
tive sign to the A 1 / 2 amplitude. It should be stressed that 
the nature of the loops, where some intermediate states 
can be put on shell, naturally leads to an imaginary part 
of the amplitude and hence one obtains complex tran- 
sition form factors. Comparison with the data implies 
a choice of phase to make the Ai/ 2 amplitude real and 
with the sign chosen in the experimental analysis. How- 
ever, once this is done, the rest of the amplitudes have 
very well determined signs and phases . In this sense we 
found that the ratio of the S1/2 to the Ai/ 2 amplitude 
was practically real and negative, and we also found that 
the ratio of the Ai/ 2 amplitude of the neutron resonance 
to that of the proton resonance was also practically real 
and of the order of -0.80, in good agreement with exper- 
iment. 

It should be noted that the signs and strengths of the 
different amplitudes are a nontrivial consequence of the 
contribution of the different channels in the photon tran- 
sition loops and of subtle interference of terms. 

Thus we can say, that the agreement with the data 
is fair when it comes to the shape of the Q 2 dependence 
and good in the ratios of amplitudes which are free of the 
global normalization. All these features together provide 
a boost to the hypothesis of the TV* (1535) as being a dy- 
namically generated resonance. This of course does not 
exclude some other components beyond those of meson 
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baryon exploited here, but the claim would be that these 
are the dominant components of the wave function and 
they show up clearly in the electromagnetic properties 
studied here. The slower experimental fall down with 
Q 2 could be an indication of the contribution of genuine 
quark components, along the lines of the work of [39| as 
we discussed in Sec. IIII Al 

The discrepancies found in the normalization of Ay 2 
for the proton deserve more attention. We have al- 
ready commented that should one use the width of the 
TV* (1535) of 90 McV of BES, or 100 MeV of MAID2007, 
the experimental values would be lowered and the agree- 
ment between theoretical results and the experiment 
would be better. In fact, the agreement of the theoreti- 
cal results for A^ 2 with the MAID2007 analysis is very 
good, as we have already noticed. But then we could 
look at the ratio R = Sf, 2 /A^, 2 and we find R = 0.6 at 

Q 2 = 0.5 GeV 2 . Experimentally, this ratio is R ~ 0.2 if 
we take an average value of Si/ 2 over the different data, 
so the discrepancies in this ratio seem to be large. Cer- 
tainly, the experimental ratio becomes much larger if we 
take the points with open triangles in Fig. [6j and then 
R ~ 0.44. This large dispersion of experimental values 
is understandable if one recalls that the contributions of 
the S1/2 term in the ep — > e'prj cross section (from where 
the data is extracted) is of the order of a few percent 
[27l ]. This, together with the experimental uncertainties 
in the normalization noted above, clearly indicate that 
large uncertainties in the experimental Si/ 2 are indeed 
present. Further improvements in Si / 2 in the future will 
reveal if the discrepancies in the ratio R pointed out here 
are deficiencies of the model or stem from present ex- 
perimental uncertainties, or both. But it is clear that 
stronger claims in favor of the theoretical model are tied 
to a better precision in this experimental ratio, thus pro- 
viding a justification for improved measurements of this 
magnitude. 

At the same time we addressed the problem of eta pho- 
toproduction on the proton and neutron with the same 
formalism. We found a cross section compatible with ex- 
periment in the ^p — > rjp reaction. This cross section 
also served to show evidence that our approach misses 
strength of the reaction at energies beyond the AT* (1535). 



This could be due to the fact that the width that we ob- 
tain for the resonance, of the order of 75 — 90 MeV, is 
smaller than the experimental one, or that the jp — > rjp 
reaction collects strength from higher energy resonances 
which are not dynamically generated and hence do not 
appear in our scheme. This issue is not settled in view 
of the large dispersion of results that one can find in the 
literature for the width of the Af*(1535), from about 90 
MeV to 350 MeV. Furthermore, the ratio of the cross 
sections of 771 — > ijn to jp — » r/p was obtained in fair 
agreement with experiment, particularly at energies close 
to the TV* (1535). 

Altogether, the information extracted in this paper 
provides support for the idea of the JV*(1535) resonance 
as being largely dynamically generated from the interac- 
tion of mesons and baryons, the dynamics of which seems 
to be well accounted for by chiral Lagrangians together 
with a proper coupled channels unitary treatment of the 
interaction, as provided by the chiral unitary approach. 

However, we also discussed that the recent study of 
[39| indicates the need for a genuine quark component of 
the JV*(1535), which could provide strength at large Q 2 
where our model, both in the nonrclativistic and rela- 
tivistic versions, still provides a faster fall down with Q 2 
than experiment. 
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